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Docket No.: 1454.1205 
IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 

In re the Application of: 

PCT National Phase of PCT/DE00/01754 

Georg DENK 

Serial No. Group Art Unit: To be assigned 

Confirmation No. 

Filed: Examiner: To be assigned 

For: COMPUTER-ASSISTED METHOD FOR THE PARALLEL CALCULATION OF THE 
OPERATING POINT OF ELECTRIC CIRCUITS 

PRELIMINARY AMENDMENT 

Assistant Commissioner for Patents 
Washington, D.C. 20231 

Sir: 

Before examination of the above-identified application, please amend the application as 

follows: 

IN THE ABSTRACT: 

Please REPLACE the Abstract originally filed with the enclosed Substitute Abstract 
attached hereto. 

IN THE SPECIFICATION: 

Please REPLACE the specification originally filed with the enclosed Substitute 
Specification. 



IN THE CLAIMS: 

Please CANCEL claims 1-10. 



PCT National Phase of PCT/DE00/01754 
Please ADD new claims 1 1-24 in accordance with the following: 

1 1 . (NEW) A computer-aided method for parallel calculation of the operating point of 
electrical circuits, comprising: 

partitioning the circuit into a number of partitions in a first step; and 
using the charging method for the parallel calculation of the individual partitions, 
a dynamic element being provided at each node of the circuit. 

12. (NEW) The computer-aided method as claimed in claim 11, wherein each node of 
the circuit is connected to in each case a predetermined value having in each case a potential 
by means of in each case one capacitance so that an operating point of the modified circuit can 
be calculated. 

13. (NEW) The computer-aided method as claimed in claim 12, wherein a capacitance 
having the same value is provided at each node of a partition. 

14. (NEW) The computer-aided method as claimed in claim 12, wherein each node of a 
partition is connected to the same potential by means of a capacitance. 

15. (NEW) The computer-aided method as claimed in claim 12, wherein a capacitance 
having the same value is provided at each node of all partitions. 

16. (NEW) The computer-aided method as claimed in claim 12, wherein 

each node of all partitions is connected to the same potential by means of a 

capacitance. 

17. (NEW) The computer-aided method as claimed in claim 12, wherein the potential is 
connected to ground. 

18. (NEW) The computer-aided method as claimed in claim 12, wherein 

the operating point of the circuit is calculated in each case with a suitable step- 
by-step change in the value of (C) of the capacitance, and 

this step is repeated until the values of the capacitances are almost zero. 

19. (NEW) The computer-aided method as claimed in claim 13, wherein each node of a 
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partition is connected to the same potential by means of a capacitance. 

20. (NEW) The computer-aided method as claimed in claim 19, wherein a capacitance 
having the same value is provided at each node of all partitions. 

21. (NEW) The computer-aided method as claimed in claim 20, wherein 

each node of all partitions is connected to the same potential by means of a 

capacitance. 

22. (NEW) The computer-aided method as claimed in claim 21, wherein the potential is 
connected to ground. 

23. (NEW) The computer-aided method as claimed in claim 22, wherein 

the operating point of the circuit is calculated in each case with a suitable step- 
by-step change in the value of (C) of the capacitance, and 

this step is repeated until the values of the capacitances are almost zero. 

24. (NEW) A computer readable medium storing a program to control a computer to 
perform a method for parallel calculation of the operating point of electrical circuits, the method 
comprising: 

partitioning the circuit into a number of partitions in a first step; and 
using the charging method for the parallel calculation of the individual partitions, a 
dynamic element being provided at each node of the circuit. 
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REMARKS 



This Preliminary Amendment is submitted to improve the form of the specification as 
originally-filed. A substitute specification and marked-up copy of the original specification are 
enclosed. No new matter is added to these documents. 

It is respectfully requested that this Preliminary Amendment be entered in the above- 
referenced application. 

If any further fees are required in connection with the filing of this Preliminary 
Amendment, please charge same to our Deposit Account No. 19-3935. 



Respectfully submitted, 



STAAS & HALSEY LLP 




700 Eleventh Street, NW, Suite 500 
Washington, D.C. 20001 
(202) 434-1500 
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SUBSTITUTE ABSTRACT 

A computer-aided method for parallel calculation of the operating point of electrical 
circuits has the circuit partitioned into a number of partitions in a first step, in which the charging 
method is used for the parallel calculation of the individual partitions. 
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SUBSTITUTE SPECIFICATION 

TITLE OF THE INVENTION 

COMPUTER-AIDED METHOD FOR PARALLEL CALCULATION OF THE OPERATING POINT 

OF ELECTRICAL CIRCUITS 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] This application is based on and hereby claims priority to German Application No. 199 
27 301.4 filed on June 15, 1999 in Germany, and PCT Application No. PCT/DE00/0 17754 filed 
on May 30, 2000, the contents of which are hereby incorporated by reference. 

BACKGROUND OF THE INVENTION 

[0002] The computer-aided simulation of electrical circuits has attained increasing importance 
in the development of very large circuits, that is to say circuits having a very large number of 
elements. It is particularly in the development of computer chips having a multiplicity, for 
example several hundred thousand transistors, that serial processing for the determination of 
the circuit quantities by a computer has been found to be unusable because of the excessive 
time consumption. 

[0003] In WO 98/24039, therefore, it is proposed to partition a large circuit and to have the 
partitions processed by different computers in each case. 

[0004] In the calculation, the operating point, that is to say the potentials of all nodes is usually 
determined first as basis for further analyses such as, for example, transient or alternating- 
current analyses. 

[0005] For the parallel calculation, an implementation of the Newton method is proposed in U. 
Wever, Q. Zheng et al.: "Domain Decomposition Methods for Circuit Simulation" (Proceedings 
of the 8th Workshop on Parallel and Distributed Simulation, PADS '94 Edinburgh, Scotland, UK, 
pp. 183-186, July 1994) and in U. Wever, Q. Zheng: "Parallel Transient Analysis for Circuit 
Simulation" (Proceedings of the 29th Annual Hawaii International Conference on System 
Sciences, pp. 442-447, 1996). The disadvantage is that convergence can only be achieved 
here when sufficiently good estimates of the operating point are available, due to poor 
convergence characteristics. As a rule, however, such good estimates can be achieved with 
difficulty or not at all in the case of large circuits. 

SUMMARY OF THE INVENTION 



Docket No.: 1454.1205 
Inventor: Georg DENK 

[0006] One aspect of the present invention is, therefore, based on the object of creating a 
computer-aided method for parallel calculation of the operating point of electrical circuits which 
ensures simple, reliable and fast calculation of the operating point of the circuit. 

[0007] In the method, the charging method, which has very good convergence characteristics, 
is used for the parallel calculation of the operating point. 

[0008] As described, for example, in H. Spiro: "Simulation of integrated circuits" (2nd edition, R. 
Oldenbourg Verlag Munich, Vienna 1990), the charging method uses the dynamic elements 
such as capacitances and/or inductances existing in an electrical circuit for calculating the 
operating point via a pseudo-transient analysis. In this process, the dynamic elements are 
charged up step-by-step to a value which, in principle, can have any magnitude, particularly to 
the value "1", starting from a value zero. 

[0009] In an embodiment according to the invention, dynamic elements such as, for example, 
capacitances or inductances having predetermined known values can be additionally installed at 
suitable places, for example at or between a number of nodes. This makes it possible to avoid 
problems in the calculation due to the number of dynamic elements being too small. 

[0010] For this purpose, the dynamic elements additionally inserted into the circuit can be 
reduced step-by-step to zero from an initial value, for example tending to infinity, for which the 
operating point calculation is trivial, so that the original circuit is again simulated. 

[0011] In a preferred embodiment of the invention, at least one node, that is to say a junction of 
at least two current paths, of the circuit is connected by a capacitance to a predetermined value 
having a predetermined potential. 

[0012] In a further development of the invention, however, a capacitance, the second terminal 
of which is in each case connected to a predetermined potential, for example to ground, can 
also be connected to each node via all partitions. This procedure has the advantage that the 
calculation of the operating point, i.e. of the respective potentials of the individual nodes, is 
trivial for the circuit for an initial value for the capacitances which tend toward infinity, due to the 
equation, which is then explicit, being solved. Changing the value for the capacitances step-by- 
step then makes it possible to change the circuit simulation by suitable new selection of the 
value for the capacitances until a calculation of the operating point of the circuit is obtained for a 
value of the capacitances tending toward zero or almost toward zero. 
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[0013] By also implementing this charging method for parallel calculation of an electrical circuit, 
a very large circuit having a multiplicity of transistors can be advantageously calculated in a fast 
and simple manner even with a smaller number of dynamic elements. Various procedures are 
conceivable for the respective new determination of the value for the capacitance with the aim 
of allowing this to go toward zero, the decision criterion being the degree of difficulty of 
calculating the operating point of the respective preceding step. 

[0014] The operating point obtained by solving a nonlinear equation for a particular value of the 
capacitances in each case can be solved, for example, iteratively by the Newton method. For 
choosing the next value for "C", the number of iteration steps necessary for the solution can 
then be used until the value drops below a predetermined value for "C". 

□ BRIEF DESCRIPTION OF THE DRAWINGS 

'ass? 

Q [0015] These and other objects and advantages of the present invention will become more 

,2 apparent and more readily appreciated from the following description of the preferred 

Nj embodiments, taken in conjunction with the accompanying drawings of which: 

* [0016] In the text which follows, the invention will be explained by means of an exemplary 
Q embodiment shown in the drawing. In the drawing, the single figure shows a flowchart which 
ilJ represents the individual method steps of the method. 

□ DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

iU 

[0017] Reference will now be made in detail to the preferred embodiments of the present 
invention, examples of which are illustrated in the accompanying drawings, wherein like 
reference numerals refer to like elements throughout. 

[0018] According to the drawing, a simulation of an electrical circuit is the starting point 1, for 
example in the circuit description language SPICE. This is partitioned in a first processing 
process 2, as described, for example, in WO 98/24039 as "Clustering method", so that individual 
partitions or parts of the circuits are obtained which can be calculated with the same degree of 
difficulty, if possible. 

[0019] Naturally, other applicable variants of partitioning are conceivable such as, for example, 
the "ratio cut method" described in N. Frohlich, B. Riess, U. Wever, Q. Zheng: "A new 
approach for parallel simulation of VLSI circuits on a transistor level" (IEEE Transaction on 
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Circuits and Systems - I: Fundamental Theory and Applications, Vol. 45, No. 6, June 1998, 
pages 601 to 613), or even an arbitrary "manual" division. 

[0020] In a further process step 3, a grounded capacitance is added to each node of the circuit, 
that is to say to a junction of at least two conductors or current paths, respectively. Naturally, it 
is also conceivable to add a capacitance, the second terminal of which is connected to a 
predetermined potential, in which case both potentials and values for the capacitances of each 
node can differ. For reasons of efficiency, stability and economy of calculation, each 
capacitance is connected to the same potential with its second terminal, for example ground, an 
identically high value CO also being selected for all capacitances in step 4 for the same reasons. 

[0021] For this value CO, the operating point is then calculated for each partition or part-circuit 
in a further step 6, in which the required coupling values, that is to say the values for the 
coupling points or interfaces of adjacent partitions being exchanged and included in the 
calculation of the operating points of adjacent partitions. 

[0022] In this process, a partition which will be called "master 1 ' in the text which follows can 
advantageously take over the control of the charging process for reasons of efficiency. The 
master then determines the initial value CO for the capacitances, and it is also conceivable to 
predetermine the initial value externally, for example by the user. This value CO is then 
transferred to all other partitions, called slaves in the text which follows. Following this, the 
operating point is calculated both in the master partition and in all slave partitions, during which 
process, naturally, the full source vector, which represents the energy sources existing in the 
circuit, is present at the circuit. 

[0023] Starting from an overall problem to be solved 

x = (m, s x , • . • , s p ) 



where m is the unknowns of the master, 

p is the number of partitions, 

Sj is the unknowns of partition i, 
the charging process leads to the system of differential equations, 

f(x,t) + D^L = 0 
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where 

D = diag (C, ... C, 0, 0) and 
t is the time. 

[0024] This system of differential equations can be solved, for example, by the implicit Euler 
method which leads to the non-linear equation 

f(x k+1 ,t k+1 ) + iD(x k+I -x k ) = 0 

where x k is the solution at time t k for k = 0, 1 , etc and h is the step t k+1 - t k . This non- 
linear equation can be solved, for example, iteratively by the Newton method 




for n = 0, 1 , etc. 

[0025] In the parallel calculation, the master then only calculates the unknowns "m", the 
unknowns "s" representing fixed areas for the master. To be able to calculate a solution of this, 
the slaves or slave partitions calculate the corrected values sf +1 in each iteration step of the 

master and report these to the master. For this purpose, the slaves must solve the system of 
non-linear equations. 

f i (sr,mr,t k+1 )+^D i ((sr i ,mr i )-(sf,m k ))=0 

where fj and Dj are the corresponding equations and matrices for the slave "i". This non- 
linear equation can also be calculated by the Newton method and not all iteration steps have to 
be necessarily performed to convergence. 

[0026] The results of the slaves are then inserted into the system equations of the master 
whereupon the master can calculate the expression mll\ . 

[0027] In a next step 7, a new value "C new " for "C" is determined by the master and the slaves 
can also make suggestions. To achieve the aim, namely a value for "C" tending to zero or, 
respectively, less than a predetermined value e, and thus a calculation of the original circuit, the 
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choice of the new value "C new " for "C" is made dependent on the difficulty of calculating the 
preceding step with the value "C 0 id" for "C". 

[0028] In this context, various procedures are conceivable, for example the choice of "C new " in 
dependence on an analysis of the number of iteration steps which the master needed for 
calculating the non-linear system of equations for the preceding value "C 0 id" for "C": 



•'old/ 



y old 



is n<n 



isn 1 <n <n 2 



C -2 is n 9 >n 

old 1 



where "n" is the number of iterations of the master and "n/ "n 2 " are the parameters 
predetermined by the user. Naturally, instead of halving or doubling "C 0 id'\ respectively, other 
strategies for reducing or enlarging "C 0 | d " are also possible. Additionally, the number of iteration 
steps needed by the slaves or slave partitions for solving their non-linear system of equations 
can also be taken into consideration, for example through the choice of 



/ 11 u 

n = max n,^n i) slavel,... ? ^n i? slavep 



i=l 



where n^ siavej is the number of iteration steps of slave j, during the i-th iteration of the 
master. It must be noted in this context that the choice of a value for "C" only affects the 
efficiency of the method not the operating point itself. 

[0029] After a value for "C" of less than or equal to a predetermined value e has been reached, 
the calculation is terminated at branch 5 and the value for C can be set to "0" in a last step. The 
operating point of the original circuit is thus found as result 8 and can be output via output units 
such as for example, a screen, printer or the like and/or stored in a memory as basis for further 
analyses of the circuit. 

[0030] Due to the method described above, an operating point of a very large electrical circuit 
can be advantageously calculated in parallel by a multiplicity of computers or processors and 
the disadvantage of known parallel types of calculation, namely lack of convergence with 
unfavorable starting values can be avoided. Due to the parallel calculation, the expenditure for 
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calculation with regard to iteration steps and setting a value for "C" step by step, can be kept 
within tolerable limits due to the distribution to a number of processors. 

[0031] The operating point thus calculated is then used as basis for further analyses, for 
example the alternating-current analysis of a circuit. 

[0032] The invention has been described in detail with particular reference to preferred 
embodiments thereof and examples, but it will be understood that variations and modifications 
can be effected within the spirit and scope of the invention. 
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MARKED-UP OF THE TRANSLATED INTERNATIONAL APPLICATION 

rDescriptionl TITLE OF THE INVENTION 

COMPUTER-AIDED METHOD FOR PARALLEL CALCULATION OF THE OPERATING POINT 

OF ELECTRICAL CIRCUITS 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] This application is based on and hereby claims priority to German Application No. 199 
27 301.4 filed on June 15, 1999 in Germany, and PCT Application No. PCT/DEOO/0 17754 filed 
on May 30, 2000, the contents of which are hereby incorporated by reference. 

BACKGROUND OF THE INVENTION 

[0002] The computer-aided simulation of electrical circuits has attained increasing importance 
in the development of very large circuits, that is to say circuits having a very large number of 
elements. It is particularly in the development of computer chips having a multiplicity, for 
example several hundred thousand transistors, that serial processing for the determination of 
the circuit quantities by a computer has been found to be unusable because of the excessive 
time consumption. 

[0003] In WO 98/24039, therefore, it is proposed to partition a large circuit and to have the 
partitions processed by different computers in each case. 

[0004] In the calculation, the operating point, that is to say the potentials of all nodes is usually 
determined first as basis for further analyses such as, for example, transient or alternating- 
current analyses. 

[0005] For the parallel calculation, an implementation of the Newton method is proposed in 
U. Wever, Q. Zheng et al.: "Domain Decomposition Methods for Circuit Simulation" 
(Proceedings of the 8th Workshop on Parallel and Distributed Simulation, PADS '94 Edinburgh, 
Scotland, UK, pp. 183-186, July 1994) and in U. Wever, Q. Zheng: "Parallel Transient Analysis 
for Circuit Simulation" (Proceedings of the 29th Annual Hawaii International Conference on 
System Sciences, pp. 442-447, 1996). The disadvantage is that convergence can only be 
achieved here when sufficiently good estimates of the operating point are available, due to poor 
convergence characteristics. As a rule, however, such good estimates can be achieved with 
difficulty or not at all in the case of large circuits. 

SUMMARY OF THE INVENTION 
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[0006] fThel One aspect of the present invention is, therefore, based on the object of creating a 
computer-aided method for parallel calculation of the operating point of electrical circuits which 
ensures simple, reliable and fast calculation of the operating point of the circuit. [According to 
the invention, this object is achieved by the method as claimed in claim 1.] 

[0007] In the method, the charging method, which has very good convergence characteristics, 
is used for the parallel calculation of the operating point. 

[0008] As described, for example, in H. Spiro: "Simulation of integrated circuits" (2nd edition, R. 
Oldenbourg Verlag Munich, Vienna 1990), the charging method uses the dynamic elements 
such as capacitances and/or inductances existing in an electrical circuit for calculating the 
operating point via a pseudo-transient analysis. In this process, the dynamic elements are 
charged up step-by-step to a value which, in principle, can have any magnitude, particularly to 
the value "1", starting from a value zero. 

[0009] In an embodiment according to the invention, dynamic elements such as, for example, 
capacitances or inductances having predetermined known values can be additionally installed at 
suitable places, for example at or between a number of nodes. This makes it possible to avoid 
problems in the calculation due to the number of dynamic elements being too small. 

[0010] For this purpose, the dynamic elements additionally inserted into the circuit can be 
reduced step-by-step to zero from an initial value, for example tending to infinity, for which the 
operating point calculation is trivial, so that the original circuit is again simulated. 

[0011] In a preferred embodiment of the invention, at least one node, that is to say a junction of 
at least two current paths, of the circuit is connected by [means of] a capacitance to a 
predetermined value having a predetermined potential. 

[0012] In a further development of the invention, however, a capacitance, the second terminal 
of which is in each case connected to a predetermined potential, for example to ground, can 
also be connected to each node via all partitions. This procedure has the advantage that the 
calculation of the operating point, i.e. of the respective potentials of the individual nodes, is 
trivial for the circuit for an initial value for the capacitances which tend toward infinity, due to the 
equation, which is then explicit, being solved. Changing the value for the capacitances step-by- 
step then makes it possible to change the circuit simulation by suitable new selection of the 
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value for the capacitances until a calculation of the operating point of the circuit is obtained for a 
value of the capacitances tending toward zero or almost toward zero. 

[0013] By also implementing this charging method for parallel calculation of an electrical circuit, 
a very large circuit having a multiplicity of transistors can be advantageously calculated in a fast 
and simple manner even with a smaller number of dynamic elements. Various procedures are 
conceivable for the respective new determination of the value for the capacitance with the aim 
of allowing this to go toward zero, the decision criterion being the degree of difficulty of 
calculating the operating point of the respective preceding step. 

[0014] The operating point obtained by solving a nonlinear equation for a particular value of the 
capacitances in each case can be solved, for example, iteratively by [means of] the Newton 
method. For choosing the next value for "C", the number of iteration steps necessary for the 
solution can then be used until the value drops below a predetermined value for "C". [Further 
advantageous embodiments of the invention are obtained from the subclaims.] 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] These and other objects and advantages of the present invention will become more 
apparent and more readily appreciated from the following description of the preferred 
embodiments, taken in conjunction with the accompanying drawings of which: 

[0016] In the text which follows, the invention will be explained by means of an exemplary 
embodiment shown in the drawing. In the drawing, the single figure shows a flowchart which 
represents the individual method steps of the method. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

[0017] Reference will now be made in detail to the preferred embodiments of the present 
invention, examples of which are illustrated in the accompanying drawings, wherein like 
reference numerals refer to like elements throughout. 

[0018] According to the drawing, a simulation of an electrical circuit is the starting point 1, for 
example in the circuit description language SPICE. This is partitioned in a first processing 
process 2, as described, for example, in WO 98/24039 as "Clustering method", so that individual 
partitions or parts of the circuits are obtained which can be calculated with the same degree of 
difficulty, if possible. 
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[0019] Naturally, other applicable variants of partitioning are conceivable such as, for example, 
the "ratio cut method" described in N. Frohlich, B. Riess, U. Wever, Q. Zheng: "A new approach 
for parallel simulation of VLSI circuits on a transistor level" (IEEE Transaction on Circuits and 
Systems - I: Fundamental Theory and Applications, Vol. 45, No. 6, June 1998, pages 601 to 
613), or even an arbitrary "manual" division. 

[0020] In a further process step 3, a grounded capacitance is added to each node of the circuit, 
that is to say to a junction of at least two conductors or current paths, respectively. Naturally, it is 
also conceivable to add a capacitance, the second terminal of which is connected to a 
predetermined potential, in which case both potentials and values for the capacitances of each 
node can differ. For reasons of efficiency, stability and economy of calculation, each capacitance 
is connected to the same potential with its second terminal, for example ground, an identically 
high value CO also being selected for all capacitances in step 4 for the same reasons. 

[0021] For this value CO, the operating point is then calculated for each partition or part-circuit 
in a further step 6, in which the required coupling values, that is to say the values for the 
coupling points or interfaces of adjacent partitions being exchanged and included in the 
calculation of the operating points of adjacent partitions. 

[0022] In this process, a partition which will be called "master" in the text which follows can 
advantageously take over the control of the charging process for reasons of efficiency. The 
master then determines the initial value CO for the capacitances, and it is also conceivable to 
predetermine the initial value externally, for example by the user. This value CO is then 
transferred to all other partitions, called slaves in the text which follows. Following this, the 
operating point is calculated both in the master partition and in all slave partitions, during which 
process, naturally, the full source vector, which represents the energy sources existing in the 
circuit, is present at the circuit. 

[0023] Starting from an overall problem to be solved 

x = (m, Sj, • . . , s p ) 



where m is the unknowns of the master, 

p is the number of partitions, 



Docket No.: 1454.1205 
Inventor: Georg DENK 



Sj is the unknowns of partition i, 
the charging process leads to the system of differential equations, 

f(x,t) + D^ = 0 



where 

D = diag (C, ... C, 0, 0) and 
t is the time. 

[0024] This system of differential equations can be solved, for example, by [means of] the 
implicit Euler method which leads to the non-linear equation 



where x k is the solution at time t k for k = 0, 1 , etc and h is the step t k+1 - 1\ This non- 
linear equation can be solved, for example, iteratively by [means of] the Newton method 



for n = 0, 1 , etc. 

[0025] In the parallel calculation, the master then only calculates the unknowns "m", the 
unknowns "s" representing fixed areas for the master. To be able to calculate a solution of this, 
the slaves or slave partitions calculate the corrected values s k+1 in each iteration step of the 

master and report these to the master. For this purpose, the slaves must solve the system of 
non-linear equations. 



where fj and Dj are the corresponding equations and matrices for the slave T. This non- 
linear equation can also be calculated by [means of] the Newton method and not all iteration 
steps have to be necessarily performed to convergence. 

[0026] The results of the slaves are then inserted into the system equations of the master 
whereupon the master can calculate the expression m*** . 



f(x k+l ? t k+1 ) + ^D(x k+1 -x k ) = 0 
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[0027] In a next step 7, a new value "C ne w" for "C" is determined by the master and the slaves 
can also make suggestions. To achieve the aim, namely a value for "C" tending to zero or, 
respectively, less than a predetermined value e, and thus a calculation of the original circuit, the 
choice of the new value "C ne w" for "C" is made dependent on the difficulty of calculating the 
preceding step with the value "C 0 i d " for "C". 

[0028] In this context, various procedures are conceivable, for example the choice of "Cnew" in 
dependence on an analysis of the number of iteration steps which the master needed for 
calculating the non-linear system of equations for the preceding value "C 0 id" for "C": 



r = 

new 



"old/ 



y old 



is n < ii 



is iij < n < n 2 



C «2 is n ? > n 

old * 



where "n" is the number of iterations of the master and "nr "n 2 " are the parameters 
predetermined by the user. Naturally, instead of halving or doubling "C 0 id", respectively, other 
strategies for reducing or enlarging "C 0 | d " are also possible. Additionally, the number of iteration 
steps needed by the slaves or slave partitions for solving their non-linear system of equations 
can also be taken into consideration, for example through the choice of 



f 



n = max 



n,J]n i ,slavel v ..,2n i ,slavep 

V i=l 



i=l 



where n i( siavej is the number of iteration steps of slave j, during the i-th iteration of the 
master. It must be noted in this context that the choice of a value for "C" only affects the 
efficiency of the method not the operating point itself. 

[0029] After a value for "C" of less than or equal to a predetermined value e has been reached, 
the calculation is terminated at branch 5 and the value for C can be set to "0" in a last step. The 
operating point of the original circuit is thus found as result 8 and can be output via output units 
such as for example, a screen, printer or the like and/or stored in a memory as basis for further 
analyses of the circuit. 
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[0030] Due to the method [according to the inventionl described above , an operating point of a 
very large electrical circuit can be advantageously calculated in parallel by a multiplicity of 
computers or processors and the disadvantage of known parallel types of calculation, namely 
lack of convergence with unfavorable starting values can be avoided. Due to the parallel 
calculation, the expenditure for calculation with regard to iteration steps and setting a value for 
"C" step by step, can be kept within tolerable limits due to the distribution to a number of 
processors. 

[0031] The operating point thus calculated is then used as basis for further analyses, for 
example the alternating-current analysis of a circuit. 

[0032] The invention has been described in detail with particular reference to preferred 
embodiments thereof and examples, but it will be understood that variations and modifications 
can be effected within the spirit and scope of the invention. 
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Description 

Compu ter-aided method for par allel calculation^ of the^operating 
point of electric al circuits 
5 ~ * 

The computer-aided simulation of electrical circuits has attained 
increasing importance in the development of very large circuits, 
that is to say circuits having a very large number of elements. It 
is particularly in the development of computer chips having a 
10 multiplicity , for example several hundred thousand transistors, 
that serial processing for the determination of the circuit 
quantities by a computer has been found to be unusable because of 
the excessive time consumption. 

15 In WO 98/24039, therefore, it is proposed to partition a large 
circuit and to have the partitions processed by different 
computers in each case. 

In the calculation, the operating point, that is to say the 
20 potentials of all nodes is usually determined first as basis for 
further analyses such as, for example, transient or alternating- 
current analyses. 

For the parallel calculation, an implementation of the Newton 
25 method is proposed in U. Wever, Q. Zheng et al . : * Domain 
Decomposition Methods for Circuit Simulation" (Proceedings of the 
8th Workshop on Parallel and Distributed Simulation, PADS '94 
Edinburgh, Scotland, UK, pp. 183-186, July 1994) and in U. Wever, 
Q. Zheng: w Parallel Transient Analysis for Circuit Simulation" 
30 (Proceedings of the 29th Annual Hawaii International Conference on 
System Sciences, pp. 442-447, 1996) . The disadvantage is that 
convergence can only be achieved here when sufficiently good 
estimates of the operating point are available, due to poor 
convergence characteristics. As a rule, however, such good 
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estimates can be achieved with difficulty or not at all in the 
case of large circuits. 
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The present invention is, therefore, based on the object of 
creating a computer-aided method for parallel calculation of the 
operating point of electrical circuits which ensures simple, 
reliable and fast calculation of the operating point of the 
circuit. 

According to the invention, this object is achieved by the method 
as claimed in claim 1. 

In the method, the charging method, which has very good 
convergence characteristics, is used for the parallel calculation 
of the operating point. 

As described, for example, in H. Spiro: "Simulation of integrated 
circuits" (2nd edition, R. Oldenbourg Verlag Munich, Vienna 1990) , 
the charging method uses the dynamic elements such as capacitances 
and/or inductances existing in an electrical circuit for 
calculating the operating point via a pseudo- transient analysis. 
In this process, the dynamic elements are charged up step-by- step 
to a value which, in principle, can have any magnitude, 
particularly to the value w l", starting from a value zero. 

In an embodiment according to the invention, dynamic elements such 
as, for example, capacitances or inductances having predetermined 
known values can be additionally installed at suitable places, for 
example at or between a number of nodes. This makes it possible to 
avoid problems in the calculation due to the number of dynamic 
elements being too small. 

For this purpose, the dynamic elements additionally inserted into 
the circuit can be reduced step-by- step to zero from an initial 
value, for example tending to infinity, for which the operating 
point calculation is trivial, so that the original circuit is 
again simulated. 
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In a preferred embodiment of the invention, at least one node, 
that is to say a junction of at least two current paths, of the 
circuit is connected by means of a capacitance to a predetermined 
value having a predetermined potential. 

5 

In a further development of the invention, however, a capacitance, 
the second terminal of which is in each case connected to a 
predetermined potential, for example to ground, can also be 
connected to each node via all partitions. This procedure has the 
10 advantage that the calculation of the operating point, i.e. of the 
respective potentials of the individual nodes, is trivial for the 
circuit for an initial value for the capacitances which tend 
toward infinity, due to the equation, which is then explicit, 
M; being solved. Changing the value for the capacitances step-by- step 

15 then makes it possible to change the circuit simulation by 
jSKSj suitable new selection of the value for the capacitances until a 

; fl calculation of the operating point of the circuit is obtained for 

a value of the capacitances tending toward zero or almost toward 



zero. 



20 



H; By also implementing this charging method for parallel calculation 

i-lf of an electrical circuit, a very large circuit having a 

\J multiplicity of transistors can be advantageously calculated in a 

;!f! fast and simple manner even with a smaller number of dynamic 

25 elements. Various procedures are conceivable for the respective 
new determination of the value for the capacitance with the aim of 
allowing this to go toward zero, the decision criterion being the 
degree of difficulty of calculating the operating point of the 
respective preceding step. 



30 



The operating point obtained by solving a nonlinear equation for a 
particular value of the capacitances in each case can be solved, 
for example, iteratively by means of the Newton method. For 
choosing the next value for W C" , 
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the number of iteration steps necessary for the solution can then 
be used until the value drops below a predetermined value for "C" . 

Further advantageous embodiments of the invention are obtained 
5 from the subclaims. 

In the text which follows, the invention will be explained by 
means of an exemplary embodiment shown in the drawing. In the 
drawing, the single figure shows a flowchart which represents the 
10 individual method steps of the method. 

According to the drawing, a simulation of an electrical circuit is 
the starting point 1, for example in the circuit description 
language SPICE. This is partitioned in a first processing process 
2, as described, for example, in WO 98/24039 as "Clustering 
method", so that individual partitions or parts of the circuits 
are obtained which can be calculated with the same degree of 
difficulty, if possible. 

Naturally, other applicable variants of partitioning are 
conceivable such as, for example, the "ratio cut method" described 
in N. Frohlich, B. Riess, U. Wever, Q. Zheng: "A new approach for 
parallel simulation of VLSI circuits on a transistor level" (IEEE 
Transaction on Circuits and Systems - I: Fundamental Theory and 
Applications, Vol. 45, No. 6, June 1998, pages 601 to 613), or 
even an arbitrary "manual" division. 

In a further process step 3, a grounded capacitance is added to 
each node of the circuit, that is to say to a junction of at least 
30 two conductors or current paths, respectively. Naturally, it is 
also conceivable to add a capacitance, the second terminal of 
which is connected to a predetermined potential, in which case 
both potentials and values for the capacitances of each node 
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can differ. For reasons of efficiency, stability and economy of 
calculation, each capacitance is connected to the same potential 
with its second terminal, for example ground, an identically high 
value CO also being selected for all capacitances in step 4 for 
the same reasons. 

For this value CO, the operating point is then calculated for each 
partition or part -circuit in a further step 6, in which the 
required coupling values, that is to say the values for the 
coupling points or interfaces of adjacent partitions being 
exchanged and included in the calculation of the operating points 
of adjacent partitions. 

In this process, a partition which will be called "master" in the 
text which follows can advantageously take over the control of the 
charging process for reasons of efficiency. The master then 
determines the initial value CO for the capacitances, and it is 
also conceivable to predetermine the initial value externally, for 
example by the user. This value CO is then transferred to all 
other partitions, called slaves in the text which follows. 
Following this, the operating point is calculated both in the 
master partition and in all slave partitions, during which 
process, naturally, the full source vector, which represents the 
energy sources existing in the circuit, is present at the circuit. 

Starting from an overall problem to be solved 
x = (m, . . . , s p ) 

where m is the unknowns of the master, 
p is the number of partitions, 
Si is the unknowns of partition i, 
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f(x,t)+D^ = 0 

where 



5 D = diag (C, ... C, 0, 0) and 

t is the time. 



This system of differential equations can be solved, for example, 
by means of the implicit Euler method which leads to the non- 
10 linear equation 

f(x k + i t k + i )+ i D(x k + i_ x k )=0 

where x k is the solution at time t k for k = 0, 1, etc and h is the 
15 step t k+1 - t k . This non-linear equation can be solved, for example, 
iteratively by means of the Newton method 

«xS* -[f x (xrSt k i+iDp^f(xr , ,t k+1 )+iD(x n k+I -x k j 

20 for n = 0, 1, etc. 



Q In the parallel calculation, the master then only calculates the 

\U unknowns w m" , the unknowns w Si" representing fixed areas for the 

master. To be able to calculate a solution of this, the slaves or 

k+l 

25 slave partitions calculate the corrected values S { in each 
iteration step of the master and report these to the master. For 
this purpose, the slaves must solve the system of non- linear 
equations . 
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r,mr i ,t k+I )+^D i ((sr i ,mr i )-(sNm k ))=0 
n 



gat 
'ft* 



ilj 



where fi and D A are the corresponding equations and matrices for 
the slave m i" . This non-linear equation can also be calculated by 
means of the Newton method and not all iteration steps have to be 
necessarily performed to convergence. 

The results of the slaves are then inserted into the system 
equations of the master whereupon the master can calculate the 
expression tn n+1 . 

In a next step 7, a new value w C new " for U C" is determined' by the 
master and the slaves can also make suggestions. To achieve the 
aim, namely a value for U C" tending to zero or, respectively, less 
than a predetermined value e, and thus a calculation of the 
original circuit, the choice of the new value w C new " for tt C" is 
made dependent on the difficulty of calculating the preceding step 
with the value w C 0 i d " for W C" . 

In this context, various procedures are conceivable, for example 
the choice of w C ne w" in dependence on an analysis of the number of 
iteration steps which the master needed for calculating the non- 
linear system of equations for the preceding value w C 0 i d " for tt C" : 
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10 



15 



where w n" is the number of iterations of the master and w ni* w n 2 " 
are the parameters predetermined by the user. Naturally, instead 
of halving or doubling w C old " , respectively, other strategies for 
reducing or enlarging u C 0 id" are also possible. Additionally, the 
number of iteration steps needed by the slaves or slave partitions 
for solving their non-linear system of equations can also be taken 
into consideration, for example through the choice of 

n > X 11 i ' slaVel > ■ ' • » X 11 i 9 slaVe P 

i=l t=l J 

where n ± , siavej is the number of iteration steps of slave j , during 
the i-th iteration of the master. It must be noted in this context 
that the choice of a value for U C" only affects the efficiency of 
the method not the operating point itself. 



After a value for *C* of less than or equal to a predetermined 
value e has been reached, the calculation is terminated at branch 
5 and the value for C can be set to "0" in a last step. The 
operating point of the original circuit is thus found as result 8 
20 and can be output via output units such as for example, a screen, 
printer or the like and/or stored in a memory as basis for further 
analyses of the circuit. 

Due to the method according to the invention, an operating point 
25 of a very large electrical circuit can be advantageously 
calculated in parallel by a multiplicity of computers or 
processors and the disadvantage of known parallel types of 
calculation, namely lack of convergence with unfavorable starting 
values can be avoided. Due to the parallel calculation, the 
3 0 expenditure for calculation with regard to iteration steps and 



GR 99 P 2008 



- 9 - 



setting a value for U C" step by step, can be kept within tolerable 
limits due to the distribution to a number of processors. 

The operating point thus calculated is then used as basis for 
further analyses, for example the alternating- current analysis of 
a circuit. 
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Patent claims 



10 



A computer-aided method for parallel calculation of the 
operating point of electrical circuits, 

in which the circuit is partitioned into a number of 

partitions in a first step, 
characterized in that 

the charging method is used for the parallel calculation 

of the individual partitions. 

The computer-aided method as claimed in claim 1, 
characterized in that a dynamic element (C, L) is provided 
and at least one node of the circuit. 



15 3. The computer-aided method as claimed in claim 2, 
characterized in that a dynamic element (C, L) is provided 
at each node of the circuit. 

4. The computer-aided method as claimed in claim 2 or 3, 
20 characterized in that each node of the circuit is connected 

by means of in each case one capacitance to in each case a 
predetermined value having in each case a potential so that 
an operating point of the modified circuit can be 
calculated. 



25 



The computer-aided method as claimed in claim 4, 
characterized in that a capacitance having the same value 
(CO) is provided at each node of a partition. 



30 6. 



The computer-aided method as claimed in claim 4 or 5, 
characterized in that each node of a partition is connected 
to the same potential by means of a capacitance. 
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7. The computer-aided method as claimed in claim 4, 
characterized in that a capacitance having the same value 
(CO) is provided at each node of all partitions. 

5 8. The computer-aided method as claimed in claim 4 or 7, 
characterized in that each node of all partitions is 
connected to the same potential by means of a capacitance. 

9. The computer-aided method as claimed in one of claims 4 to 
10 7, characterized in that the potential is connected to 

ground . 

10. The computer-aided method as claimed in one of claims 4 to 
8, characterized in that 

15 - the operating point of the circuit is calculated in each 

case with a suitable step-by- step change in the value of 
(C) of the capacitance, and 

this step is repeated until the values of the 
capacitances are almost zero. 

20 

11. A computer program product which can be loaded into a main 
memory of a computer system, with a software code for 
carrying out the method according to one of the preceding 
claims when the computer program product is running on a 

25 computer system. 



12. 



A data carrier with a computer program product as 
claimed in claim 11. 



GR 99 P 2008 
Abstract 

Computer-aided method for parallel calculation of the operating 
point of electrical circuits 

The invention relates to a computer-aided method for parallel 
calculation of the operating point of electrical circuits, in 
which the circuit is partitioned into a number of partitions in a 
first step, in which the charging method is used for the parallel 
calculation of the individual partitions. 



The single figure is the main drawing. 
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Patent claims 

1. A computer-aided method for parallel calculation of the 
operating point of electrical circuits, 

5 - in which the circuit is partitioned into a number of 

partitions in a first step, 
characterized in that 

the charging method is used for the parallel calculation 
of the individual partitions, a dynamic element (C, L) 
10 being provided at each node of the circuit. 

2. The computer-aided method as claimed in claim 1, 
characterized in that each node of the circuit is connected 
to in each case a predetermined value having in each case a 

15 potential by means of in each case one capacitance so that 

an operating point of the modified circuit can be 
calculated. 

3. The computer-aided method as claimed in claim 2, 
characterized in that a capacitance having the same value 
(CO) is provided at each node of a partition. 



4. The computer-aided method as claimed in claim 2 or 3 , 
characterized in that each node of a partition is connected 

25 to the same potential by means of a capacitance. 

5. The computer-aided method as claimed in claim 2, 
characterized in that a capacitance having the same value 
(CO) is provided at each node of all partitions. 

30 

6. The computer-aided method as claimed in claim 2 or 
5, characterized in that 
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each node of al-1 partitions is connected to the same 
potential by means of a capacitance. 

7. The computer-aided method as claimed in one of claims 2 to 

5, characterized in that the potential is connected to 
ground . 

8. The computer-aided method as claimed in one of claims 2 to 

6, characterized in that 

the operating point of the circuit is calculated in each 
case with a suitable step-by- step change in the value of 
(C) of the capacitance, and 

this step is repeated until the values of the 
capacitances are almost zero. 

9. A computer program product which can be loaded into a main 
memory of a computer system, with a software code for 
carrying out the method according to one of the preceding 
claims when the computer program product is running on a 
computer system. 

10. A data carrier with a computer program product as claimed in 
claim 9. 
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Ich beanspruche hiermit auslandische Prioritatsvorteile 
gemass Abschnitt 35 der Zivilprozessordnung der 
Vereinigten Staaten, Paragraph 119 aller unten ange- 
gebenen Auslandsanmeldungen fur ein Patent oder 
eine Erfindersurkunde, und habe auch aile Auslands- 
anmeldungen fur ein Patent oder eine Erfindersurkun- 
de nachstehend gekennzeichnet, die ein Anmelde- 
datum haben, das vor dem Anmeldedatum der 
Anmeldung liegt, fur die Prioritat beansprucht wird. 



f hereby claim foreign priority benefits under Title 35, 
United States Code, §119 of any foreign application(s) 
for patent or inventor's certificate listed below and have 
also identified below any foreign application for patent 
or inventor's certificate having a filing date before that 
of the application on which priority is claimed: 
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Prior foreign appplications 
Prioritat beansprucht 



19927301.4 DE 



(Number) 
(Nummer) 



(Number) 
(Nummer) 



(Number) 
(Nummer) 



(Country) 
(Land) 



(Country) 
(Land) 



(Country) 
(Land) 



15.06.1999 

(Day Month Year Filed) 

(Tag Monat Jahr eingereicht) 



(Day Month Year Filed) 
(Tag Monat Jahr eingereicht) 



(Day Month Year Filed) 
(Tag Monat Jahr eingereicht) 



Priority Claimed 



^ □ 
Yes No 
Ja Nein 



□ 

Yes 
Ja 



□ 

Yes 
Ja 



□ 
No 
Nein 



□ 
No 
Nein 



Ich beanspruche hiermit gemass Absatz 35 der Zivil- 
prozessordnung der Vereinigten Staaten, Paragraph 
120, den Vorzug aller unten aufgefuhrten Anmel- 
dungen und falls der Gegenstand aus jedem Anspruch 
dieser Anmeldung nicht in einer fruheren 
amerikanischen Patentanmeldung laut dem ersten 
Paragraphen des Absatzes 35 der Zivilprozefcordnung 
der Vereinigten Staaten, Paragraph 122 offenbart ist, 
erkenne ich gemass Absatz 37, Bundesgesetzbuch, 
Paragraph 1.56(a) meine Pflicht zur Offenbarung von 
Informationen an, die zwischen dem Anmeldedatum 
der fruheren Anmeldung und dem nationalen oder PCT 
internationalen Anmeldedatum dieser Anmeldung 
bekannt geworden sind. 



I hereby claim the benefit under Title 35. United States 
Code. §120 of any United States application(s) listed 
below and, insofar as the subject matter of each of the 
claims of this application is not disclosed in the prior 
United States application in the manner provided by 
the first paragraph of Title 35, United States Code, 
§122, I acknowledge the duty to disclose material 
information as defined in Title 37, Code of Federal 
Regulations, §1. 56(a) which occured between the filing 
date of the prior application and the national or PCT 
international filing date of this application. 



PCT/DE00/01754 



(Application Serial No.) 
(Anmeldeseriennummer) 



(Application Serial No.) 
(Anmeldeseriennummer) 



30.05.2000 

(Filing Date D, M, Y) 
(Anmeldedatum T, M, J) 



(Filing Date D.M.Y) 
(Anmeldedatum T, M; J) 



anhangig 

(Status) 

(patentiert, anhangig, 
aufgegeben) 



(Status) 

(patentiert, anhangig, 
aufgeben) 



pending 

(Status) 

(patented, pending, 
abandoned) 



(Status) 

(patented, pending, 
abandoned) 



Ich erklare hiermit, dass alle von mir in der vorliegen- 
den Erklarung gemachten Angaben nach meinem 
besten Wissen und Gewissen der vollen Wahrheit 
entsprechen, und dass ich diese eidesstattliche Erkla- 
rung in Kenntnis dessen abgebe, dass wissentlich und 
vorsatzlich falsche Angaben gemass Paragraph 1001, 
Absatz 18 der Zivilprozessordnung der Vereinigten 
Staaten von Amerika mit Geldstrafe belegt und/oder 
Gefangnis bestraft werden koennen, und dass derartig 
wissentlich und vorsatzlich falsche Angaben die Gul- 
tigkeit der vorliegenden Patentanmeldung oder eines 
darauf erteilten Patentes gefahrden konnen. 



I hereby declare that all statements made herein of my 
own knowledge are true and that all statements made 
on information and belief are believed to be true, and 
further that these statements were made with the 
knowledge that willful false statements and the like so 
made are punishable by fine or imprisonment, or both 
under Section 1001 of Title 18 of the United States 
Code and that such willful false statements may 
jeopardize the validity of the application or any patent 
issued thereon. 
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VERTRETUNGSVOLLMACHT: Ais benannter Erfinder 
beauftrage ich hiermit den nachstehend benannten 
Patentanwalt (oder die nachstehend benannten 
Patentanwalte) und/oder Patent-Agenten mit der 
Verfolgung der vorliegendert Patentanmeidung sowie 
mit der Abwicklung aller damit verbundenen Geschafte 
vor dem Patent- und Warenzeichenamt: (Name und 
Registrationsnummer anfuhren) 



POWER OF ATTORNEY: As a named inventor, I 
hereby appoint the following attorney(s) and/or 
agent(s) to prosecute this application and transact all 
.business in the Patent and Trademark Office 
connected therewith, (fist name and registration 
number) 



Customer No. 21171 



And I hereby appoint 



Telefongesprache bitte richten an: 
(Name und Telefonnummer) 



Direct Telephone Calls to: (name and telephone 
number) 

Ext. . 



Postanschrift: Send Correspondence to: 

Staas & Halsey LLP 
700 Eleventh Street NW, Suite 500 20001 Washington, DC 
Telephone: (001) 202 434 1500^an^Fac§imile (001) 202 434 1501 

Custom 



H 





VolleJ^aTrle/^ej^einzigen oder ursprunglichen Erfinders: 

Dr. GEORG DENK 



Full name of sole or first inventor: 

Dr. GEORG DENK 



Unterschrift des Erfin3ers^ 

1 fo3c 



Datum 



Inventor's signature 



Date 



k XlM 



&.Q<Df 



Wohn: 



ROSENHEIM, DEUTSCHLAND 

Staatsangehorigkeit 



4tt 



Residette 

ROSENHEIM, GERMANY 



DE 



Citizenship 

DE 



Postanschrift 

FARRENPOINTSTR. 5B 



Post Office Addess 

FARRENPOINTSTR. 5B 



83026 ROSENHEIM 



83026 ROSENHEIM 



Voller Name des zweiten Miterfinders (falls zutreffend): 



Full name of second joint inventor, if any: 



Unterschrift des Erfinders 



Datum 



Second Inventor's signature 



Date 



Wohnsitz 



Residence 



Staatsangehorigkeit 



Citizenship 



Postanschrift 



Post Office Address 



(Bitte entsprechende Informationen und Unterschriften im (Supply similar information and signature for third and 
Faile von dritten und weiteren Miterfindern angeben). subsequent joint inventors). 
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